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One and two photon wave functions are obtained by projection onto a basis of simultaneous 
eigenvectors of the position and number operators. 



The current interest in entanglement and its applica- 
tion to quantum communications has rekindled the de- 
bate on the nature of the photon wave function 
The photodetection amplitude has been identified with 
the real space photon wave function in the discussion of 
down conversion experiments 0] , a choice that can be jus- 
tified by it's relationship to photon counting which can 
localize the photon |a, tD| • In the standard formulation of 
quantum mechanics the real space wave function is the 
projection of the state vector onto an orthonormal basis 
of eigenvectors of a Hermitian position operator. How- 
ever, it has been claimed since the early days of quantum 
mechanics that there is no position operator that defines 
such a basis for the photon. Here we will briefly review 
our recent work on the construction of a photon position 
operator and obtain a photon wave function by project- 
ing onto its eigenvectors. 

Attempts to arrive at a photon position operator and 
its associated basis of localized states go back to the early 
days of quantum mechanics. Pauli stated that the nonex- 
istence of a density for the photon corresponds to the fact 
that the position of a photon cannot be associated with 
any operator in the usual sense . Based on definitions 
of center of mass, Pryce found the k-space photon po- 
sition operator rp = iV — ik/2fc^ + k x S/fc^ where Sj 
are the 3x3 spin 1 matrices, k is a wave vector, and 
This operator does not have commuting 
components, and thus three spatial coordinates cannot si- 
multaneously have a definite value. In 1949 Newton and 
Wigner sought rotationally invariant localized states and 
the corresponding position operators. They were success- 
ful in the case of massive particles and zero mass particles 
with spin and 1/2, but found for photons "no localized 
states in the above sense exist" . This result is widely 
quoted as a proof of the nonexistence of a photon posi- 
tion operator. It has been proved that there is no pho- 
ton position operator with commuting components that 
transforms as a vector • 

Recently we have constructed a position operator with 
commuting components that is not rotationally invariant 
[l^ , does not transform as a vector , and thus is con- 
sistent with the previous work. Description of a localized 
state requires a sum over all k and a localized photon can 
have definite spin in the k-direction, that is it can have 
definite helicity, but it cannot have definite spin along 
any fixed axis. It is the total angular momentum (AM) 
that has a definite value along some specified direction in 
space 0, 0| . The position eigenvectors are not spheri- 



cally symmetric, instead the y h ave a vortex structure as 
is observed for twisted light [la| ■ Compared to the New- 
ton Wigner position operators for which transformation 
of a particle's spin and position are separable, the photon 
position operator must incorporate an additional unitary 
transfomation that reorients this vortex. 

Maxwell's equations are analogous to the Dirac equa- 
tion when written in terms of the Riemann-Silberstein 
field vector F = E±zcB [H ^ where E and B are 
the electric and magnetic fields. This suggests that the 
photon is an elementary particle like any other, and 
that Maxwell's equations provide a first quantized de- 
scription of the photon. The use of the positive fre- 
quency Riemann-Silberstein vector as a photon wave 
function has been thoroughly studied [l3 . l20j . If a 
field 'S'^^/^) such as F that goes as fc^/^ is used as 
wave function, a metric factor is required in the 
scalar product. The real space squared norm then goes 
as Jd^rJd^r'^W^)* (j.) . *(i/2) (r') / |r - and thus 

its integrand cannot be interpreted as a local number 
density [T^. Since the photon has no mass, it has been 
suggested that there is no photon number density, only 
energy density 20]. However, the Landau-Peierls (LP) 
wave function, whose absolute value squared has 

been interpreted as photon number density was investi- 
gated as early as 1930 0,113. It has the disadvantage 
that its relationship to electric current density and the 
electromagnetic fields is nonlocal in real space |l7ll23[l24l |. 
It is possible to define a biorthonormal basis with a local 
scalar product that involves the eigenvectors of an opera- 
tor and its adjoint jl^. This formalism has recently been 
applied to pseudo-Hermitian Hamiltonians that possess 
real spectra jl^. We will show here that such a basis 
provides an interesting alternative to explicit inclusion 
of a metric operator when working with electromagnetic 
fields. 

In this letter one and two photon wave functions will 
be obtained by projection onto a basis of position opera- 
tor eigenkets. Our work on the photon position operator 
will first be reviewed and the properties of biorthonor- 
mal bases will be outlined. The position eigenkets will 
be obtained in the Heisenberg picture (HP). We will then 
derive photon wave functions from quantum electrody- 
namics (QED) by projecting the state vector onto simul- 
taneous eigenkets of the photon position operator and the 
number operator. Finally we will discuss the relationship 
of these projections to the photodetection amplitude and 
other real space descriptions of the photon state in the 



2 



recent literature. 

We start with a discussion of the photon position 
operator. A k-space position operator with com- 
muting components and transverse eigenvectors in the 
spherical polar 9 and directions was introduced in 
Ref. It was generalized in Ref. to 

allow for rotation about k through the Euler an- 
gle x{0,<i>) to give r("'X) = D (k'^iVk-"') D-^ where 
D — exp ^— iS-kx^ cxp {—iSj,4>) exp {^1320) ■ The uni- 
tary transformation D rotates k from the z-axis to an 
orientation described by the angles and (f>, while the 
transverse vectors x and y are rotated first to and 
and then about k through x to give the unit vectors 



exp (— iAx) ( + iX(f> 



with helicity A = ±1. The similarity transformation k" 
results in eigenkets proportional to k" where we are in- 
terested in a = and ±1/2 as discussed above. The 
k-space position operator, 



k k X S k- S 
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fc^ fc^ /c^ 



(^cot0-Vx), (1) 



has transverse 3-vector eigenkets satisfying 



?(a,x)^("). 



(k) = nV^, 
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for a photon with helictiy Ai at ri. In Eq. ||2Jl the func- 
tional dependence distinguishes k-space from r-space, 
while subscripts denote eigenvalues and Cartesian com- 
ponents. There is a remarkable analogy between the last 
term in this k-space position operator and the r-space 
vector potential of a magnetic monpole where the Eu- 
ler angle x corres pon ds to a change of gauge. This was 
explored in Ref. [13|. It turns out that the analogy is 
primarily mathematical, and the last term in Eq.l^ does 
change the physics. The spin and orbital AM of a pho- 
ton are not separable Q . However, the z-component of 
the total AM operator commutes with the position op- 
erator and and this allows r'^^'X)^ ^j^e helicity operator 
k • S, Jz = {—ip X V -I- S) to have simultaneous eigen- 
vectors with eigenvalues ri, Ai, and hj^ for integral j^- 
The the first three terms of r^^/^''^-' are the Pryce position 
operator, rp, whose components do not commute. The 
last term gives the position operator commuting compo- 
nents, dictates that r'"'^) transform as a vector only for 
rotations about the z-axis as can be seen from Eq.(67) 
of Ref. and fixes jz for a given r'^^'X). The quantum 
numbers {ri, Ai} index the basis states for a given j^- 

A biorthonormal basis of one photon position eigen- 
kets will now be obtained. For a = the operator r(O'X) 
is self adjoint, has real eigenvalues, and defines a sin- 
gle orthonormal basis as is usual in quantum mechanics. 
For fields, a = 1/2 and the position operator is not self- 
adjoint, rather it is pseudo-Hermitian. The biorthonor- 
mal pairs, {ipm (t>n} , of eigenkets of a pseudo-Hermitian 



operator and its adjoint satisfy |25ll2 

d\ij„)^OnHn), Ot |0„) = O; |0„) , 
= rjO-q^^, {4>n\lkrn) = ^n,m, 



(3) 



where 77 is a metric operator, li p = ^ is the pos- 
itive square root of 77, then o = pOp^^ is Hermitian. 
To apply this formalism to the photon we take -q = k 
and a = —1/2. Then o = fC'-x) is Hermitian and the 
eigenvectors of O = r(-i/2,x) ^oid — r'^^/^'^-' form a 
biorthogonal pair that go as 1/Vk and Vfc as required 
by QED for the vector potential and the electromagnetic 
fields respectively. Eqs. (|3J| then give Eq. and 

j(-i/2.x)t ^ A;f(-i/2.x)y^-i ^f(i/2,x) 



(a) 
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where 5'^ is the 3-dimensional Dirac (5-function and we can 
interchange a with —a. By multiplying Eq. (0) forf^'^'-^-' 
by p^^ where p = fc^^^ and inserting p^^ p^^ — 1 between 
f and ■0 to obtain the r(='=i/2'A:) eigenvector equations it 
can be proved that the real eigenvectors, ri, are preserved 
by the similarity transfomation to the biorthogonal basis. 

The time dependence is determined by the Hamilto- 
nian H + Hq with H = ^ hkca\^ a'^^.a where the zero 

point terms Hq — hkc/2 which are unaffected by 

the photon state will be omitted here. The operator 
ak,A annihilates a photon with wave vector k and he- 
licity A. The operators and their eigenkets are time de- 
pendent in the HP |27l |. Using the unitary time evo- 
lution operator U (t) — exp (^—iHtj , the position op- 
erator, given by Eq. ^ in the Schrodinger picture, be- 
comes f^]?' = C/t (t)f("'A)[/ (i) in the HP with eigenkets 
W (t) |ri, Ai) given by 

^r?k,j C^' ^) = ^"ek,Ai,j exp (-zk • n + ikct) /W (5) 

in the k-space representation. Equivalently we can de- 
scribe the 1-photon position eigenkets by defining the 
operators 

V-r" Ai,i (t) = E ^"sk.Aij-akai exp (-ik • ri + ikct) /W 

k 

(6) 



and the kets 



(7) 



where |0) is the vacuum state. The field operators are 
E — —dA/dt and B = V x A where the vector potential 
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operator in the Coulomb gauge can be written as Ar {t) 



;(+) 



(t) + a[ ^ /V2 with 



A(+) (t)^CY^ fc^^/^ekAflkA exp (ik-r- ikct) /Vv 



(-) 
r,A 



it) 



:(+)t 

r,A 



{t) , C = yJhjctQ, and cq 



where A. 

the permittivity and c the speed of Ught in vacuum. 
The 1-photon operators given by Eq.® are simply re- 
lated to the vector potential and electric field operators 



through (k,t) = A\-XJC and V^^^^^. = 

A general state vector in which the number of photons 
and their wave vectors are uncertain can be expanded as 



:(l/2) 



I*) ^co|0>+^Ck,A<AlO) (8) 

k,A 

+ ^ XI v/-^k,A;k',A'Ck,A;k',A'ak,A4',A' |0) + - 
■ k,A;k',A' 

where cq = (0|^') , Ck,A = (0|ak,A|*), Ck,A;k',A' = 

Ck',A';k,A = (0 lok.Aftk'.A'l ) and A/^c^AikVA' = 1 + 

(5k,k"^A.A'- Division by 2! corrects for identical states ob- 
tained when the {k, A} subscripts are permuted while 
\fM /2 normalizes doubly occupied states. The one 'pho- 
ton real space wave function, equal to the projection of 
this state vector onto an eigcnket of f ^j?' is 

(r, t)^j2 (^S J*) = E ^k,A^^l;- (k, t) (9) 

A k,A 

= E ''k,Aek,A,j^" "^xp (ik • r — ikct) /Vv 

k,A 

where we have used Eqs. JT)), © and ||SJ). If a = this 
is the LP wave function, '4''^*'^ (r, t). The vector potential 

A(+) (r,i) =:C*(-i/2) (j.^^) 
determines the positive frequency fields 

(r, t) = -^A(+) (r, t) = icC*(i/2) (j.^ ^) ^ 

(r,t) = V X A(+) (r,i), 
which satisfy Maxwell's equations. The photodetection 



wave function is E^+' (r, t) 
scalar product 



E(+) (r,t) 



. The 



(^1^) = E/ (r,t) 



k,A 



has a local integrand and |ci| is the probability for 1- 
photon in state 



In k-space the 1-photon LP and field wave functions 
identically predict probability |ck,A|^ to measure mo- 
mentum Kk. In real space the LP wave function leads 
to a positive definite density j^''-''-' (r, i)| . When us- 
ing fields is replaced with the biorthonormal pair 
1^(1/2)^ ^(-1/2) I States with definite photon energy or 
angular momentum can have a definite k jl4l |. and the 
relationship between their description in the LP and the 
biorthogonal bases is trivial, even in real space. How- 
ever, for position eigenkets and real space wave functions 
in general these two bases are not so simply related. Ac- 
cording to the competeness relation in Eq.Q a = 1/2 
and —1/2 can be interchanged and the two options aver- 
aged to give the real density 

n (r, t) = Re |*(i/2)* (j., t) ■ *(-i/2) (j., t)} 
= Re {ieoE(-) (r, t) ■ A(+' (r, t) /h'j 

which is local but not positive definite, and thus it is 
not a true probability density. This can be seen from 
the following example: If \'^) is a 1-photon state that 
includes only wave vectors ki and k2 both with helicity 
A where Cki.A = Cka.A = 1/V2 then 




The cosine term can exceed the spatially uniform time 
independent terms due to the Vk factors, leading to 
negative values. It gives zero if an integral over all 
space or all time is performed, explaining why the 
scalar product is unaffected by the similarity transfor- 
mation. In an experiment that integrates over a long 
enough time or a large enough spatial volume, use of 
and the biorthonormal pair ^(='=1/2) ig equivalent. 
The density ieoE'^^ • A'^+^/2?i -I- c.c has appeared be- 
fore. The classical linear and angular field momenta 
can be written as Ylj j d^'^'i^^E-jOAj /h 1^. This form 
can be applied to optical beam AM calculations within 
the paraxial approximation |28l| . The number operator 

^ — *£oX]j-^j a[^^ /2h + h.c. transforms as the ze- 
roth component of a four-vector and satisfies a continu- 
ity equation [i^ . The one photon density n (r, t) equals 
(^I^ \n\ and its integral over all space is time indepen- 
dent consistent with the pair j'J''-'^/^), ^f^"^/^^} forming 
a basis as implied by Eqs. I@J. Action, which is of this 
form, has appeared in calculations of laser linewidth (s^ ■ 

For two photons we can project jvP) onto the 2-photon 
real space basis i>'f^l^ (t) - {t') |0) . Use of Eq.© and 



Ok.A.flk'.A' 



= (5k.k'(5. 



k.k'OA.A' 



then gives the correlation am- 
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plitude 



2\V 

k,A;k',A' 

* * ik-r—ikctik'-r'—ik'ct' 



+ek',A',i6k,Aj 



gZk-r' — zfcci' ^ik' -r—ik' ct 



(10) 



which becomes a two photon wave function if we set 
t' — t. It is a symmetric linear combination of products 
of one photon wave functions in agreement with Refs. 

and 0. The one and two photon amphtudes given 
by Eqs. ® and ((Tni) are consistent with the use of the 
detection amphtude as a wave function 0, S IE IE ■ 
The " two-photon quantum state in coordinate space" ob- 
tained by taking the Fourier transform of the k-space 
probabihty amphtude in Ref. |3l| is an example of re- 
cent use of the a — LP wave function. In cither case, 
the sum over all n-photon terms for all positions and 
helicities provides a complete basis and thus "encodes 
the maximum total knowledge describing the system" 
as required by Schrodinger and discussed in Ref. 
For example, either form can be used to transform from 
the photon number basis to the quadrature basis if all 
nonzero n-photon amplitudes are known. 

In summary, we have reviewed our previous work 
where it is demonstrated that a photon position oper- 
ator does indeed exist. Because photon spin and or- 
bital angular momentum are inseparable, its eigenvec- 
tors have a vortex structure like twisted light. We pro- 
jected the QED state vectors onto simultaneous eigenvec- 
tors of this position operator and the number operator 
in two different ways: If all k's are weighted equally the 
Landau-Peierls wave function is obtained. This gives a 
positive definite probability density, but a nonlocal rela- 
tionship to fields and current sources. A biorthonormal 
basis consisting of eigenkets proportional to the vector 
potential and electric field results in a real local density, 
ieoE^~) • A*^+)/2?i-f c.c, which is not positive definite. 
Both of these wave functions have played a role in recent 
analyses of two photon entanglement P, 0, |^ ■ The two 
photon wave function is a symmetrized product of one 
hoton wave functions in agreement with Refs. and 
When all photon numbers are allowed for, either ba- 
sis provides a complete description of the quantum state 
of the electromagnetic field, equivalent to the QED state 
vector. 
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